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Abstract: We derive finite time, dynamical radiation reaction effects from lightfront QED. 

Combining the Hamiltonian formalism with a null background field, the calculation is per- 

^H formed with no approximation other than the usual coupling expansion of perturbative 

"^ QED. The literature contains many proposals for the correct classical equation describing 

^ a radiating particle; we take the classical limit of our results and identify which equations 

are consistent with QED. 
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1 Introduction 

The motion of a particle is influenced not only by external forces but also by the particle's 
own emission of radiation. Momentum conservation implies that the particle recoils when 
it radiates, an effect called 'radiation reaction' (RR). 

A rough order-of-magnitude estimate, equating the electron rest energy to the work 
done by an electric field over the classical electron radius, suggests that RR becomes signif- 
icant at electric field strengths of roughly ^ lO^^V/m (though see below). This is two orders 
of magnitude higher than the Sauter-Schwinger limit at which nonperturbative QED effects 
come into play [1, 2]. Despite this, and the implied difficulty in observation, understanding 
RR presents one of the oldest and most frequently revisited problems in electrodynamics. 
Let us recall why. 

Beginning with the coupled classical equations for motion for a particle, electromagnetic 
fields and external forces, it is possible to integrate out the field variables and write down an 
equation for the orbit of a radiating particle. This is the well known Lorentz-Abraham-Dirac 
('LAD') equation [3-5], which has two unusual features. The first is that a divergence arises 
in its derivation, but this can be removed by renormalising the mass. The second feature 
of LAD is that it is third order in derivatives and admits unphysical runaway solutions in 
which even a free particle can spontaneously accelerate to the speed of light. Due to this. 



a great deal of work has over the years gone into deriving the 'correct' equation of motion 
for a classical, radiating particle which avoids the problems of LAD. 

Runaways are nonperturbative in the electromagnetic coupling, being given by 1/e^ 
terms [6]. The simplest way to avoid them is therefore, given a particular system, to 
solve LAD perturbatively. To any given order the solutions are free from runaways (and 
pre-acceleration). To obtain an equation which describes such physical solutions, one can 
reduce the order of LAD from third to second by recursively substituting the equation into 
itself, which eliminates higher derivative terms. To obtain a tractable equation one must 
additionally truncate this expansion to some order in the coupling. The Landau-Lifshitz 
equation (LL) [7] is the first order truncation, and is probably the most commonly employed 
equation to describe classical RR. The results of [8] support LL as the effective equation 
describing physical, i.e. non-runaway solutions to first order in the coupling^. Many further 
proposals for classical equations exist, and it is fair to say that there is no absolute consensus 
in the literature over which should be used to describe classical physics [9-11]. 

Perturbative quantum electrodynamics (QED) is the best tested theory of fundamental 
physics [12]. Starting in QED, we should be able to derive RR effects, and then take the 
classical limit to provide some definite answers. At first glance it seems straightforward 
to extend the above description of RR to quantum field theory. One begins with a state 
describing an electron, evolves this in time through, say, an external field which excites 
the state, causing photon emission, and then measures the momentum of the electron as a 
function of time. This should show effects due to the electron recoiling as it emits. Funda- 
mentally, though, quantum field theory is a multi-particle theory; an accelerated electron 
can radiate photons which can produce electron-positron pairs, and therefore particle num- 
ber is not conserved. The quantum notion of 'an electron's recoil' then becomes ambiguous 
since there is a varying number of electrons in the system. 

A formal way around this problem is to restrict to the regime in which there is only a 
single electron. In perturbation theory (in particular, in the Furry picture, see below), if 
one begins with a single electron state, then there remains a single electron in the system 
to order a — e^/(47r/i); zeroth order terms contain only acceleration effects due to external 
forces, while order a terms describe photon emission and self-energy, both of which are 
reminiscent of classical RR [6, 13]. At order o:^, two-photon emission contributes, but so 
does, e.g., pair production [14-16]. To extend the idea of quantum RR to higher orders, 
one can work in a parameter regime for which pair production has much lower probability 
than multi-photon emission, see [17, 18] and below. 

In a previous paper, we explained in general which diagrams contribute to lowest order 
RR effects in the S^-matrix of QED [13]; these are indeed photon emission and the electron 
self energy. It is only the inclusive combination of the two which yields a physical, measur- 
able, IR-finite observable; this will be reviewed below. We showed explicitly, for a certain 
class of background fields, that one recovers known asymptotic results for classical RR from 
our QED results, in the limit /i ^ 0: asymptotically, and to lowest order, all the classical 



^This result is sometimes misquoted as being exact, but it is clearly stated in [8] that higher orders are 
dropped. Further, LAD and LL differ already at order e^, see below. 



equations we will consider agree, and are consistent with Larmor's formula (for quantum 
corrections to which see [19, 20]). 

The purpose of this paper is to derive finite time RR effects from QED, which re- 
quires going beyond S^-matrix elements and instead investigating the dynamics of states 
at finite time. This is most natural in a Hamiltonian formalism, which we combine here 
with lightfront field theory [21, 22], in which quantisation is performed on null hyperplanes 
and x+ = x^ + x^ is the time co-ordinate. Going from S^-matrix elements to finite time 
dynamics brings its own challenges; renormalisation in the Hamiltonian formalism, which 
we encounter below, is made difficult due to a lack of explicit covariance, a problem com- 
pounded on the lightfront due to the appearance of nonlocal, momentum-dependent coun- 
terterms [23]. Thus the investigation presented here goes somewhat beyond our previous, 
S^-matrix based, calculation [13]. 

We will not try to derive a general equation from QED, but will instead consider a 
particular choice of background field. The main reason for this is that we wish to present 
a calculation with the smallest possible number of approximations and assumptions. The 
topic of RR is old, the literature vast, and there are many different approaches in play, 
see [24-26] for reviews and references. In our calculation, we will use the coupling expan- 
sion of QED, but all other parameters will be treated exactly, thus eliminating potential 
ambiguities. We work to first nontrivial order in the coupling, and explain the extension to 
higher orders. Our interest here is in theory, the question of when higher order terms are 
phenomenologically important will be addressed elsewhere. The approach we use is to com- 
pute averages of appropriate position/momentum operators, before taking their classical 
limits in order to compare with the predictions of various classical equations. 

Our choice of background is a null field, or plane wave, of arbitrary temporal profile. As 
we will see, this background is particularly amenable to a lightfront treatment and allows 
us to easily distinguish between different classical equations. It is also the background 
lying behind much of the work on 'strong field QED', which studies the use of intense 
laser light, or strong magnetic fields, for investigating physics both within the standard 
model and beyond. See [24, 27, 28] for recent reviews of this active topic. The intense 
laser fields which are now, or soon will be, available present a huge potential for observing 
RR [24, 29, 30]; it has been estimated that RR effects are actually measurable at 10^^- 
lO^^W/cm^ [31-33], corresponding to field strengths which are significantly lower than in 
the coarse estimate above, and which will be reached by next generation laser facilities. 

This paper is organised as follows. In Sect. 2 we review various proposals for classical 
equations describing a radiating particle, and summarise their predictions for motion in 
a plane wave background. In Sect. 3 we perform the lightfront quantisation of QED in 
this background, and then construct and renormalise the electron momentum operator. 
The calculation of quantum and classical RR effects in the momentum and position of a 
particle, at finite time, is presented and discussed in Sect. 4. Our conclusions are presented 
in Sect. 5. We take c = eg = 1 throughout, but h^ 1 unless otherwise stated. 
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Table 1. The radiation reaction force for a particle in an external field / := eFQ^t/m. The left 
hand column shows the antisymmetric tensor form of the force. This antisymmetry implies the 
mass-shell condition x^ = 1, using which we simplify R^ in the right hand column. 



2 Classical radiation reaction 

We begin by collecting various proposals for equations describing a classical radiating par- 
ticle in an external electromagnetic field Fext- Our interest is not in the properties of these 
individual equations per 5e, but in comparing their predictions with the classical limit of 
QED. We therefore refer the reader to the original articles, cited below, for more details. 
Most of the proposed classical equations take the form 



2 e^ 
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ext-^^^ 



(2.1) 



in which a dot denotes a derivative with respect to proper time r, and i?^ describes the 
radiation reaction force. We write / := eF^y^i/m from here on. Table 1 lists the forms 
of R^ (in part using a compact notation in which Lorentz indices should be read from 
left to right) for the equations of Lorentz- Abraham-Dirac (LAD) [3-5], Landau-Lifshitz 
(LL) [7], Eliezer/Ford-0'Connell (EFO) [34, 35], Mo-Papas (MP) [36] and Herrera (H) [37]. 
The equation originally proposed by Eliezer [34] was rederived by FO for a particle with 
structure, rather than a point particle [9, 10, 38, 39]. Some of the differences between LAD, 
MP and LL were investigated in [40]. Our final classical equation is due to Sokolov (S) [41]. 
In this case, the velocity x^ is not proportional to the momentum g^, so the equation has 
a different form than (2.1). See Appendix A. From here on we focus on the equations in 
Table 1, and state corresponding results for S. 

We will solve the above equations using a classical analogue of perturbation theory 
in the Furry picture, a technique long used in strong- field-QED [42, 43]. (The quantum 
expansion is described in Section 3.) The idea is simple: we treat radiation reaction as a 
perturbation to the Lorentz force. We rewrite (2.1) as 

o2 



x^ 



rx, + 



2 e' 



-R^ 



(2.2) 



3 Aum 

in which the coefficient of i?^ is the classical electron radius, e^ /{Aum) 2^ 3 fm. Noting that 
e now appears only in the RR term, we expand the orbit in powers of e^ (for simplicity of 
presentation; the appropriate dimensionless parameter is given below) 

x{t) = xo(t) + e^X2{T) + e^X4(T) + . . . , (2.3) 



in which subscripts indicate the power of e. We plug this into (2.1) and solve the equation 
in powers of e^. To zeroth order, the equations of motion (2.2) are 

^0 = f(xo)xo , (2.4) 

which is the Lorentz force equation for a particle, orbit xq, moving in a background field 
F^^*. (This is the exact Lorentz force equation, not that equation expanded in the coupling.) 
Higher order terms Xn in the expansion (2.3) correspond to modifications of the Lorentz 
orbit due to RR, i.e. due to the electron's recoil [13]. While this expansion allows for a 
clear separation of Lorentz force and recoil effects, it is most useful when one can solve the 
Lorentz force equation exactly, as this allows higher order terms to be calculated analytically 
as well. This will be the case for our background field, to which we now turn^. 

2.1 Classical radiation reaction in null fields 

Our chosen background is a plane wave. This is constructed from a lightlike (null) wavevec- 
tor kj^ and a transverse polarisation vector a|J^, so ka' = 0. For the null vector we have 
kjj, = cjn^, in which uj is an inverse- length scale, for example a central frequency, and we 
choose coordinates such that nx = x^ + x^ = x+, which is lightfront time [21, 22]. The field 
strength depends only on dimensionless invariant phase (j) :— kx — ujx^ as 

eF-*(,/.) = k^a'M - a'^i<P)k^ , (2.5) 

The polarisation vector and electric field are related by a'^{(j)) \— eE^{(j))/uj. Throughout, 
a dash is a derivative with respect to (j). We refer to both (j) and x+ as lightfront time. 
We consider pulses, so that E^{(j)) is either nonzero only in a finite (/)-range, or vanishes 
asymptotically, but is otherwise arbitrary. The field's spacetime structure is sketched in 
Fig. 1. All massive particles enter and leave a plane wave at the same lightfront time. Once 
they leave the pulse, they can never return to it. Further, taking (j) as the time coordinate, 
spacetime acquires a 'band structure', being cleanly separated into regions 'before', 'during' 
and 'after' the pulse. No such separation is possible when using x^ as the time coordinate, 
which manifests in properties of both the classical and quantum theories, see below. We 
now proceed to solve the classical equations of motion (2.2) in a plane wave. 

Zeroth order: Lorentz force 

To zeroth order, we need to solve (2.4), the Lorentz force equation in a plane wave. The 
solution is well known, and follows from first observing that kx^ — 0, since fcFgxt = 0, 
implying that kx^ is conserved. So, if the particle has momentum p^ when it enters the 
pulse then, with 00 = /cxq, 

Ur) = -r . (2.6) 

m 



^In our chosen background, exact solutions are actually available for LL [44] and S [41]. Since we are 
interested in comparing to perturbative QED, we need only the perturbative expansions of these results. 







x+ = 



Figure 1. A pulsed plane wave has finite extent in x+, but infinite extent in the remaining 
directions. All particles, independent of their momentum, enter and leave such a pulse at the same 
lightfront times, here x+ = and x+ = ^j , respectively. There is no such symmetry in the usual 
time parameter x^. 



The particle's proper time is, to zeroth order, proportional to lightfront time. Using this, 
(2.4) becomes a linear differential equation which can be solved immediately by exponenti- 
ation. The momentum tt := uixq of a particle moving under the Lorentz force is 



TT =p-a(0o(T)) + 



2kp 



k 



(2.7) 



in which a^ is the integral of the electric field strength^. 






(2. 



Here and below, the lower limit of the integral can be shifted to the point at which the field 
turns on; this is — oc for an asymptotically switched pulse, or can always be chosen to be 
= for a compactly supported pulse, as in Fig. 1. To accommodate both options, we will 
frequently drop the lower limits on (f) integrals, taking the relevant limit to be understood. 
From here on p^ is always an initial momentum and tt^ always depends on p^ as in (2.7). 
The solution (2.7) may be integrated directly with respect to r to find the position of the 
particle. Note that a = a((/)o(T)), i.e. a function of 00, which is a function of r. This brings 
us to an important and slightly subtle point regarding the comparison between classical 
and quantum results, which we must discuss before calculating the first RR contribution to 
the orbit. 



Reparameterisation 

The classical equations above are parameterised by proper time r. Solving them yields 
X = x(t), x(t) = drx{r) and so on. In QED, and QFT in general, no reference is made 



^No gauge potential is used in the classical calculations. 



to a particle worldline; since QFTs are multi-particle theories, a global time coordinate is 
used to parameterise state evolution [45, 46]. Hence, any dynamical quantity obtained from 
QED will be parameterised by the time coordinate which defines the quantisation surfaces^. 
In our case this will be x+ := x^ + x^, which is natural given Fig.l. We therefore need to 
parameterise our classical solutions in terms of the full x+, or, to maintain covariance the 
full (/) (and not r or ^q). This can be done exactly for, at least, LL in a plane wave [44], 
and in perturbation theory for all the considered equations. 

To zeroth order, (j) — (j){) — kpr/m, so it is simple to eliminate r in order to parameterise 
the momentum in terms of light front time (f). We have 

mx{<P) = p - a{ct>) + '^P^i'^)^-^^'^^) + o{e^) ^ 7r(</>) + 0{e^) , (2.9) 

To compare the orbit with QED we can use (2.9) to write, to zeroth order, 

^^f + Oie')^^n{<t>) + Oie'), (2.10) 

d(p (j)Q kp 

following which a ^-integration yields the orbit. At higher orders, inverting (/)(t) to obtain 
t((/)) brings in more complicated functional dependencies since <j) is not constant when RR is 
accounted for [44]. (The implied symmetry breaking of null translation invariance provides 
a potential signal for measuring RR [32].) 

Even to zeroth order, there is no closed form parameterisation of the Lorentz orbit 
using instant-form time x^. Mathematically, this is because x^{r) has no simple inverse, so 
we cannot explicitly eliminate r in favour of x^. Physically, the reason is that the chosen 
background singles out x^ as a preferred direction, not x^, as in Fig. 1. 

First order: radiation reaction 

We are now ready to solve for X2, the first correction to the Lorentz orbit due to RR. 
Inserting (2.3) into (2.2), the equation to be solved is 

X2 = f{(p0)x2 + (p2f\(po)xo + R{xo) • (2.11) 

6 4:7Tm 

One may proceed as before: solve for 02 by taking the inner product with fc^, following which 
(2.11) becomes a standard differential equation which can be solved for X2- Alternatively, 
one can begin with the {x~ ,x^} components (which are simple), and integrate directly with 
respect to r or 00, and then change variables to 0; the final component, x+ then follows 
from the mass-shell condition. Either way, the result is that the momentum of the radiating 
particle is, to order e^ and as a function of 0, 

mm = m + l~J^^ (^(-o(^)) - ^^M^^M,) + oie^) . (2.12) 

Either of the methods above can be extended directly to higher orders. It remains to insert 
the RR forces from Table (1) into (2.12). We summarise these results below. 



^It would though be interesting to reformulate our later calculations using worldline QED [47, 48]. 



2.2 Classical predictions 

Define the dimensionless parameter A by 

A:= 



2 e^ kp 



(2.13) 



3 47r 771^ 

Inserting the expressions from Table (1) into (2.12) to obtain the momentum q (equal to 
mx in all cases except S), one finds to order e^ 



q{ct>) = vr(^) + An'{ct>) + A L^ a'\^) fvr(^) - !^M!^fc 




m^ 



kp 



(2.14) 
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For the position (which follows from integration except in S, where q ^ mx) we find 
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kp^ - A.'( 

d'? LAD 



(2.15) 



The classical solutions differ in their transverse and longitudinal components. LAD, LL 
and EFO give the same result to order e^, whereas MP, H and S predict a different result, 
as they do not contain the derivative A7^^ In this context, we note that LL is usually 
obtained by performing a reduction of order on LAD; applying the same reduction to EFO 
also yields LL [26]. Reducing MP one finds H. Further, reducing MP as an equation for 
momentum g, one finds S, at least in a plane wave background. 

All the equations predict the same final momentum at = oc, which is equal to the 
momentum when the particle leaves the pulse; writing rr = 7r(oc) from here on^, this is 



q{oo) 



7r + 



A 



m' 



jd^ a'\^) (vr(^) - ^fc 



kp 



(2.16) 



This can also be derived from Larmor's formula for the total emitted radiation. The equa- 
tions also agree, at all times, on the behaviour of the momentum component fcx, which is 
conserved in the case of the Lorentz force but with lowest order RR is 







(2.17) 



^TT / p in general [49, 50], though attention is usually restricted to pulse shapes which do not give 
vacuum acceleration, in which case it — p. All our results hold irrespective of the pulse shape. 



From this non-conservation of kx, we can extract the typical scale of RR effects in a plane 
wave. If we scale out 77 = ^^ as a typical magnitude of the electric field, and imagine that 
the pulse is short so that the remaining (j) integral in (2.17) is of order unity, we find 

^ ~ 1- Ar?2 + 0(e^) = l-lf^,v' + 0{e') . (2.18) 

kp 3 4:7T m^ 

(Ary^ is an appropriate, dimensionless, expansion parameter for RR in a plane wave.) For 
a 10 PW laser system, we have 77 ^ 100 and optical frequency huj ^ 1.24 eV. Assuming a 
head-on collision between the beam and particle with gamma factor 7, we find 

lf^,v'-lO-S, (2.19) 

suggesting that RR effects reach, for example, 10% of Lorentz force effects at around 7717 = 
0.5 GeV. At 7717 ^ 5 GeV, (2.19) reaches unity, and higher order corrections must be taken 
into account (see e.g. [44] for the exact solution of LL in a plane wave.) This suggests that 
RR may indeed be measurable at the 10 PW laser facilities currently under construction. 
For more detailed estimates and predictions of RR effects we refer the reader to [31-33]. 
While higher gamma factors can reduce field strengths required for observing RR, they also 
take us toward the quantum regime. Hence it is time to consider quantum effects. 

We will compare (2.14) and (2.15) to the classical limit of QED, specifically to ex- 
pectation values of the electron momentum and position operators in the limit h ^ 0. 
Calculating these expectation values is the focus of the remainder of this paper. 

3 Lightfront quantisation of scalar QED 

Our aim is to derive classical RR from the finite time dynamics of quantum states. We 
therefore use the Hamiltonian formalism. The calculation we will perform mirrors that in 
the classical theory, where we began with an electron and solved its equations of motion. 
Here we begin with single electron state and solve the Schrodinger equation for its time 
evolution. Also as above, we treat the 'Lorentz component' exactly. For related calculations 
using the Hamiltonian approach and expectation values, see [51], 'in-in' calculations [52- 
54] and position shift calculations [55]. A comparison with S^-matrix approaches is given 
below, and position calculations in our background are presented in Sect. 4.2. 

Consider first a particle in an external plane wave, without radiation. We saw above 
that while the corresponding Lorentz equation can be solved exactly, the solution has no 
explicit parameterisation in terms of x^. Parameterisation in terms of = fcx, or lightfront 
time, is on the other hand natural and straightforward. This is reflected in the quantum 
theory of a particle in a plane wave: despite the basics of this apparently simple theory 
having long been known [56], it is only in recent years that progress has been made in the 
instant form canonical quantisation of the theory [57, 58], and even then only for a specific 
choice of plane wave. The situation is very different if one quantises on the lightfront [21, 22]. 
Quantisation then proceeds in analogy with the free theory, and has recently been used to 
clarify long-standing ambiguities regarding the 'effective electron mass' in a plane wave [59]. 



Further, when QED interactions are added, hghtfront quantisation proceeds as for ordinary 
QED on the hghtfront, as first shown in [60]. 

We consider scalar QED (sQED) from here on, as spin effects wih in any case drop 
out in the classical limit. In order to perform the calculation of interest we must set up 
the theory, regulate and renormalise. sQED comprises the gauge field A^, describing the 
photon, and a complex scalar field $, describing the scalar electron and positron. Including 
an additional background A^^*, we have the action 



f 



^x - -F^^F^^ + {D^^)^D^'^ - ^$t$ ^ (3.1) 



in which D^ = d^^-i^A^^i^A^^^ , F^^ = d^Ay — djjA^ and the background obeys Maxwell's 
equations in vacuum, d^F^^^ — 0. 

We will quantise on null hyperplanes of constant = fcx, on which the background field 
depends. This allows us to retain some explicit covariance [61]. We can, though, always 
choose our co-ordinates such that, as above, kx — uj{x^ + x^) = cjx^, the usual hghtfront 
time direction. With this choice the remaining coordinates are 'transverse' x^ — {x^,x^} 
and 'longitudinal', x~ — x^ — x^ . These will be denoted with sans-serif fonts, so x = 
{x^,x~}. Corresponding momentum components are p = {p±,P-} with p± = \{p^ ^Ps)- 
Integrals over these variables are written 

oo 

/*<./d.-dV. /.p../^/||, ,3.) 



The derivation of the hghtfront Hamiltonian follows that for sQED without background, so 
we highlight only the important steps. For a pedagogical discussion see [21]. For simplicity 
we set h — 1 and reintroduce it only before taking the classical limit. As usual, we use 
hghtfront gauge, A'^ = 2A_ — 0, for both the dynamical and background field. In hghtfront 
quantisation, only the transverse components of A^ are dynamical fields; the longitudinal 
component A^ is determined by A^ and the current J_ according to^ 

^+ ~ 2d. 2dl • ^^-^^ 

The dynamical fields are then also transverse to the background's propagation direction, 
just like the physical fields of the background itself. This is one reason for orientating our 
quantisation surfaces with the laser direction [61]. A convenient choice of lightfront-gauge 
potential for the background field is then eAJ^^{x) = ciii{(l))-> which is easily verified to 
give the correct field strength (2.5). The usefulness of this choice will soon be clear. The 
Hamiltonian can now be written down. 



H =^ fdx Aj{id^fAj + |P^$p + m^l^p , 



+ ^JdxefA^-e^A^A^\^\^ + 



. , (3.4) 



e 
J 



2'-{id- 



Vl V2 



^Inverting d- requires a prescription for dealing with zero modes [21, 22, 62], see also below. 
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in which and from here on A_^ = d^A^/2d_ [21], the background-covariant derivative is 
T^/j, — djj, ^ ia^ and the background current is j^, 

3^ = i$tp^$ _ i(P^$)t$ = ^\id^ - 2a^)$ . (3.5) 

(Note that J_ = j_.) Apart from the presence of V^ rather than 5^, (3.4) is the ordinary 
Hamiltonian of hghtfront sQED. We wih quantise in the Furry picture [63]. This is a par- 
ticular choice of interaction picture (which reduces to the usual one when the background 
vanishes) defined by the separation of the Hamiltonian into 'free' and 'interacting' parts. 
The 'free' theory is described by the first line of (3.4), Hq. It comprises free photons, and 
the scalar field interacting with the plane wave. This theory can be solved exactly [60], 
like the classical Lorentz equation can be solved exactly. The second line of (3.4) is the 
'interacting' part of the Hamiltonian, containing Vi and V2 which are respectively linear 
and quadratic in the coupling. These are the usual lightfront vertices of sQED but with the 
free matter current replaced by (4.21); the three- and four-point vertices, both including 
instantaneous-photon interactions due to A^, and finally the instantaneous-scalar interac- 
tion. This 'interacting' part of the Hamiltonian is treated in perturbation theory. In terms 
of position space Feynman diagrams, one has the usual vertices, the usual photon prop- 
agator, and a 'dressed' fermion propagator describing the scalar's propagation within the 
background. 

Operators evolve in the Furry picture under the action of the 'free' Hamiltonian H^. 
The transverse gauge field therefore has the usual mode expansion of lightfront field theory, 

Aj{x) = U\ a,(l)e-^^^ + a](l)e^^^ , (3.6) 

in which ^+ = l^/(4l_) as usual, so that t^ = 0, on shell; all particles are on-shell in lightfront 
perturbation theory [21]. The fields obey the lightfront commutation relations [22], 

[aii\),a]{\')]^2l_{2nf5^{\-\')5ij, 

with e the sign function. Defining a+(l) = /±tt±(l)/2/_, the constrained field A^ can be 
written 

A^{x) = fd\ a+(l)e-^^^ + a^l)e^^^ , (3.8) 

and one can then write down a covariant expression for the mode commutators: 

[a,(l),at(|')] = -2l_{27rf5^\ - I') (r?,. - hiji±ll^^ ^ -2l_{2n)H{\ - 1')^. (3.9) 
The scalar field operator, under the action of i7o, has the mode expansion 

$(x) = y"dp 6(p)^p(x) + dt(p)^_p(x) , (3.10) 
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in which the mode functions cpp are Volkov solutions [56], i.e. solutions to the Klein-Gordon 
equation in a background plane wave, 



(Pp{x) = exp 



ipx — I 



2pa — a 
2kp 





(3.11) 



with p^ = (p^ + 7Ti^)/(4p_) so that p^ = tti^, on-shell. The mode operators are the 
(quantised) initial data on the hyperplane = 0, when the background turns on. b^ and 
S create on-shell scalar electrons and positrons with initial momentum p^. One-particle 
states evolve under Hq to carry the momenta tt^ of the classical theory, which is seen in 
the following property of the mode functions 

^/i^p(^) = 7r^(0)(/^p(x) . (3.12) 

TT^ is of course also on-shell, 7r+ = (tt^ + 7Ti^)/(4p_), and 7v_ = p_ is conserved under Hq- 
evolution, as was the case classically when considering only the Lorentz force. The scalar 
field obeys the commutation relations 

[6(p),6t(q)] ^ [d(p),dt(q)] ^ 2p.{2nfS^p - q) . 
We now turn to the states. These evolve in the Furry picture according to 

i-^\7P;x+)F^HF{x+)\7P;x+)F, (3.14) 

in which Hp is the Furry analogue of the 'interacting Hamiltonian in the interaction picture', 

Hp = T:e'^f ^0(^1 + V2me-'y ^0 = (^i + 1^2)L=(3.6), *=(3.io) (3-15) 

i.e. one inserts the 'free field' mode expansions (3.6) and (3.10) into the second line of (3.4). 
We normal order throughout. Our interest lies in the case that the initial state describes a 
single electron. The state in the far past is therefore 

I in) = I dp g{p)b\p)\0) , Jdp |5(p)|' = 1 , (3.16) 

in which | ) is the lightfront vacuum annihilated by the a(l), 6(p) and d(p), and where ^(p) 
is a wavepacket, normalised as shown, such that |( in | in)p = 1. The state at later lightfront 
time x+ is then 



|'0;x+) = 7^ exp 



f ds HFis) 



in), (3.17) 



in which we make the usual assumption of asymptotic switching to evolve the initial state 
up to the fully interacting electron state before it enters the background field. 

This completes our quantisation of the theory. We now want to calculate the expecta- 
tion values of the electron momentum operator P^ in the evolved state (3.17), 

(pp = (V';x+|p;|V';x+), (3.18) 

and then take the classical limit in order to compare with the classical results in Sect. 2. 
First we need to write down the momentum operator. 
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3.1 The momentum operator 

The energy-momentum tensor in our theory is 

r^, = --L/^ = F^^F% + {D^^)^D,^ + {D,^)^D^^ - g^X , (3.19) 

and the totalmomentum of the system is [21, 22] 

P^ = Jdx T_^ , (3.20) 

in which P+ is the Hamihonian (3.4). We are not interested in the total momentum, but 
in that of the electron. We must therefore separate T^^^ into a piece which describes the 
electron, and a piece which describes the emitted radiation. This is trivial in the free theory, 
and also asymptotically in which the theory again becomes free^. At finite time, though, it 
is not obvious how to make such a separation, since the physical electron is a composite of 
the matter field and a cloud of photons [66-69]. A related classical problem is the subject 
of energy balance, the question of how one separates the electromagnetic fields into bound 
and radiated parts, see [70, 71] for detailed discussions. 

There are several conditions which will help us to identify the relevant part of T-^. 
First, the decomposition must be gauge invariant. Second, it must give the correct free- 
field limit. These two constraints (a third follows) suggest the natural electron-photon split 
T^^ = T^u + T;I, in which [72] 

These are the usual energy-momentum tensors for the electromagnetic field, and the min- 
imally coupled scalar. This decomposition is gauge invariant and has the correct free-field 
limit, implying we should take 

P; = jdxTl^, P; = jdxTZ,. (3.22) 

These operators are explicitly dependent on lightfront time due to the background field. To 
order e^, antiparticles do not contribute to our calculation (a common feature of lightfront 
quantisation [21, 22]), and so we drop the positron modes d{p) in all our expressions. 
Similarly, some order e^ terms in the operator drop out due to normal ordering. The 
contributing part of the electron momentum operator can then be be written as a sum of 
two terms, P^ = P/)^ ^ + P/)^ ^ ^ the first (second) being order zero (one) in the coupling: 



pW(^) = y'dp7r^(</.)fet(p)6(p), 

P(i)(,^) = eldpdl 6t(p - l)6(p) (at7r(,^)fc^ - ^^^aj, (o) 



r0 l-K 

kp — kl 



(3.23) 



''Modulo IR problems [64, 65], but IR divergences drop out of our expectation value, see [13] and below. 
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This is almost, but not quite, the operator we are looking for. Our third constraint is that 
P^ should yield the correct sQED result: an electron with momentum p^, unexposed to 
external forces, should always have momentum pj^. This means that we should find 

{P;)^P,- (3.24) 

Our first task is therefore to calculate (P^) in ordinary sQED without a background field, 
and check that we obtain the correct answer (3.24). 

3.2 On-shell renormalisation 

Regularisation and renormalisation are necessary in the following calculations. We use 
dimensional regularisation to control UV-divergences. (For classical RR in different di- 
mensions, see [73].) Following [74], the extra dimensions are placed into the transverse 
directions, so that 

fd^l^ -^ ^^' [(Th , (3.25) 

where n = 2(1 — e) and fi is the introduced mass scale. Transverse dim reg has the benefits 
of affecting neither the structure of the quantisation surfaces nor our chosen background; 
the plane wave is always homogeneous in the transverse directions, another reason for 
orientating our quantisation surfaces with the laser direction. See [75] for an application 
of transverse dim reg. Zero modes should be regulated using cutoffs, principle values or 
otherwise [21, 22, 76]. We display the regulators only when necessary, but they are in place 
throughout. In the following calculation, dim reg is sufficient to take care of divergences. 

We set a^ = to return to ordinary lightfront perturbation theory. We then begin with 
an initial electron (3.16), and calculate the expectation value of its momentum at subsequent 
times, (3.18) with P^ as in (3.23) but a^ = 0. We will present the 'in-background' version of 
this calculation in more detail later, and so we skip here to the final result. In the limit that 
the wave packet is strongly peaked around an initial momentum p^, we find to order e^. 



/De\ ,2 [a\^P~^U 2 2kplp\f Ip \ 

2e2 f 2kp-klf Ip kp \ 1/ 2 2Mp\ 

~i^J 2ip Y^^'w^'w^'j'Tpy^ ~^r) 



(3.26) 



k/j, 



The first line comes from P^^', tfie second line from P^^' . Evaluating the transverse integrals 
in dim reg, the longitudinal integrals contribute only finite terms. We find 



Pa)=P, 



27Te \ m? 



kp ^ 



(r 



27Te \m'^ 



(3.27) 



where ellipses denote finite terms. From here on we write down only divergent terms; finite 
terms are also removed by our renormalisation condition (below). The physical result p^ is 
multiplied by 'one plus a divergent quantity', and we have a second divergence proportional 
to kj^. This second divergence is typical of lightfront field theories, in that it depends on the 
momentum component kp [21-23]. The result (3.27) is though covariant, and the second 
divergence depends only on the angular variables in A:^, as expected from [61]. 
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Removal of the divergences proceeds as follows. Recall that S'-matrix elements are 
renormalised by adding counterterms to the Hamiltonian/Lagrangian. In our calculation, 
though, no order e^ term in the interaction Hamiltonian contributes to the order e^ ex- 
pectation value. Counterterms in the interaction Hamiltonian therefore cannot remove the 
divergences in this case. Instead it is the composite operator P^ (and, below, other oper- 
ators) which we need to renormalise. We will nevertheless encounter some results familiar 
from lightfront perturbation theory. 

To proceed, we need a renormalisation condition. Given that a free electron of mo- 
mentum pjj, should always have momentum j)^, none of the order a terms in (3.27) can 
be physical, and must be removed. We therefore take our renormalisation condition to be 
(3.24). It seems natural to call this a Hamiltonian analogue of on-shell renormalisation. 
The first divergence in (3.27) is proportional to p^, i.e. to the free result. It can therefore 
be removed by multiplicative renormalisation: 

1 + ^)f;. (3,28) 

In fact, the divergence comes entirely from the lowest order operator P^ , and hence only 
this part needs to be renormalised; this is equivalent to (3.28), to order e^. We turn to 
the second divergence, proportional to A:^. This can be removed by mass renormalisation. 
Using transverse dim reg, the mass shift^ in lightfront quantisation is [21, 23] 

5m = -^ . (3.29) 

ZTve 

Observe that to this order, the only affect of a mass shift on the expectation value would 
be the change 

V,^V, + 1^^ K =Pf^- -^K + 0{e^) . (3.30) 

which would precisely cancel the A:^ divergence in (3.27). Such a shift is affected by replacing 
m ^ m — 6m in the Lagrangian term of P^ (equivalent to adding a 0^ counterterm to P^), 
see (3.21). The contribution to the expectation value is precisely that required, 

m6m , 
kp 

cancelling the remaining divergence (and finite terms) in (3.27). The expectation value of 
the renormalised operator is, finally, 

(P;)renorm. ^ P^^ + C?(e^) , (3.32) 

which is the desired result^. We have identified the renormalised momentum operator which 
yields the finite, physical electron momentum to order e^. We now return to the theory 
with a background field. 



{PP ^ {Pp - -—k, , (3.31) 



^Not to be confused with the 'intensity dependent mass shift' in a plane wave, for which see [32, 59]. 

^(P^) is also divergent, and to obtain a finite result P^ can mix with P^ [77, 78], giving (P^) -^ 
(P^) — ^{P^) = + 0{e^) , the expected result. Another interpretation is that the split (3.21) is 
corrected by quantum effects, and this requires redefining the momentum operators by mixing P^ and P^: 
after transferring the p^ divergence from P^ to P^, and renormalising the electron mass, both operators 
yield finite expectation values. 
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Figure 2. The Furry- Fey nman diagram describing the zeroth order (Lorentz-force) contribution to 
the average electron momentum, (free)^. Double lines indicate the background-dressed propagator. 

4 Radiation reaction at finite time 

We now have everything required to calculate the expectation value (3.18); we have the 
state (3.17), the momentum operator (3.23), and we know how to renormalise. The above 
renormalisation yields finite results also in-background; no 'new' UV divergences arise, see 
also [79]. To order e^, the expectation value comprises the following four terms 

{PP = {m\p(%n) 

x~^ x~^ x+ y 

+ {m\ JdyVi{y)Pl,^^ ldzVi{z)\m) -2Re{m\Pl,^^ J dy J dzVi{y)Vi{z)\m) (4.1) 

x+ 

+ 2Im(in|pWydy+yi(y+)|in), 

which we write as, respectively, 

{PP = (free)^ + (emission)^ + (loop)^ + (operator^^))^ . (4.2) 

These names refer to the origin and physical meaning of each term, as we now describe. 

Free 

(Free) is the contribution from the 'free' theory without QED interactions. It corresponds 
to scattering- without-emission, e~ ^ e~, in which the electron is accelerated by the back- 
ground, but the electron does not emit. It should therefore yield the Lorentz force compo- 
nent. From (3.16) and (3.23), we find 

(free)^ = (in|Pi°)|in) = ^dp \g{p)\\^{<t>) ^ 7r^(0) , (4.3) 

where in the final step we have taken the limit in which the wavepacket is sharply peaked; 
this is indeed the Lorentz force result. See Fig. 2 for the associated diagram. 

Loop and emission 

(Loop) is the tree-level/one-loop cross term from scattering-without-emission, i.e. a con- 
tribution from the self energy of the electron. (Emission) corresponds to single photon 
emission from the electron, which in a plane wave background is called nonlinear Compton 
scattering [42, 43, 80-85]. See the diagrams in Fig. 3. To calculate the emission and loop 
terms, we need the result, writing dK = dl 9{kp — kl), 

jdy+Vr{y+)\m)^ejd<Prj^^^ 5(pKat(l)6t(p _ |)| o)e' / ^ . (4.4) 
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Pf^ 



■■P'^ X Pm 



pm 



■p^^ 



-p^ 



+ 



Pf^ 



p. 



Pa X Pa 



■■P^ 



M 



Figure 3. Furry- Fey nman diagrams for two of the O(e^) contributions to the average electron 
momentum, and to radiation reaction. In the first line, the cross term from scattering-without- 
emission, and in the second line, nonlinear Compton scattering, at tree level, mod squared. 



We find that the two terms differ only in their vector structure, and a relative minus sign. 
(This is a direct consequence of the optical theorem.) (Loop)^ contains the Lorentz force 
(no recoil) result 7r^(0) while (emission) contains p'^ the electron's momentum following 
photon emission in nonlinear Compton [13, 80-85]: 



Pu 



TTn 



-lu + 



l<<t^) 



kp 



ki^^ 



(4.5) 



Reinstating /i, we find 



(loop + emission)^ = -^^ / dl' [7r^((/)) - p^{(l))\ j 



d0id(/)2 7^ — ^e 
kp — kl 



h J kp — kl 

<t>i 



(4.6) 



The reason we group the loop and emission terms together is that while both are individually 
infra-red (IR) divergent [49], their sum is IR finite [13]. Ths is because, unlike individual 
^-matrix elements, our expectation value is an inclusive observable (no assumption is made 
about the 'final' state) [86-88]. The potentially IR-divergent term in (loop) is removed by 
(emission), due to the vector structure tt^ — P^- See [89] for related comments regarding 
the Schwinger mechanism. 

The sum (loop + emission) does though contains a UV divergence. We can separate 
out the divergent term using integration-by-parts to rewrite the ^-integrals as follows 



(loop + emission)^ 
,kp — kl 



H- 



hkp 



[7r{<t>)-p'{4>)]^ 



7r{<P) - 'SI 



kl' 



lirict^) 



d</.i 



llTl 



(4.7) 



kpixf 



exp 



In 



h J kp — kl 



The boundary terms cannot be dropped: they lead both to UV divergences which must be 
renormalised and to finite, physical terms which must be retained. They are also essential 
for distinguishing between different equations in the classical limit. (Boundary terms were 
also crucial for the comparison between LAD and QED in [55].) Expanding the modulus 
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we obtain three terms, 
(loop + emi.™.)„ = e^-JM''^ I'^*** " '''(*'l(^^^l(^)' <"■**' 

The first hne, (4.8), comes from the 'boundary-boundary' terms of the two integrals, and 
is UV divergent. Regularising, the momentum integrals in (4.8) can be performed exactly, 
and one finds 

2 
Q/ Q/ 777, 

(loop+emission)^ = -^^^M) + ^^^/^ + ^^-^^ + ^^'^O) ' (^'^l) 

The two divergences have the same form as those without background. The first is propor- 
tional to the 'free'-theory momentum 7r^(0), the second is proportional to A:^ and depends 
on kn = kp. The same multiplicative and mass renormalisations as above remove (4.8) and 
with it the two divergences. The remaining two terms, (4.9) and (4.10) are UV finite, and 
unaffected by the renormalisation. 

Operator 

Finally, (operator^-^^), is the contribution from the interacting part of the electron momen- 
tum operator, P^ K This contribution vanishes in the asymptotic limit (see below) and 
does not have a standard Feynman diagram description. Again using (4.4) we find 






2 r AU r I 

-{2kp - kl)L,^i,7Ti - ki^7T{(f))L7Tl 



_ i_ r Itt 
h J kp — kl 

01 



(4.12) 



The above expressions (minus the divergent term (4.8)) for the free, loop, emission and 
operator terms are lowest order quantum recoil effects. See [90] for numerical investigations. 

4.1 The classical limit 

The simplest way to take the classical limit follows from noting that photon momentum 
has no classical analogue, whereas wavenumber does: we therefore write /^ = hk'^ with fc|^ 
the outgoing photon's wavenumber. h can then be taken to zero in, (4.9), (4.10) and (4.12), 
and the integrals performed to obtain the classical result. 

To better understand the physics of this limit, though, it is helpful to first highlight the 
dependence of the quantum results on the longitudinal momentum component kl = hkk\ 
Longitudinal momentum is conserved by the Lorentz force, but this symmetry is broken 
by photon emission [32]; this non-conservation is seen explicitly when one goes from the 
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solution of Lorentz to the solution of e.g. LAD or LL, see above and [44]. Define the ratio 
of outgoing momenta by ht := hkk' /kp' ^ which may be equivalently written 

ht = = . (4.13) 

kp — hkk' kp — kl 

We then change integration variables from /_ to t, and we scale out the longitudinal mo- 
mentum from l^ by defining r^ — -^lii-, which amounts to a second change of variables, 
l^ ^ r^ — {kp/kl)l^. We illustrate this using the transverse component of (4.10), which 
may be equivalently written 

OO ^2 

in which the symmetry of the integrals allows us to replace exp with cos. Factors of h are 
removed from the exponents and now occur only in the combination 1 + ht. The classical 
limit therefore corresponds to Taylor expanding these factors, assuming ht <^ 1. From 
(4.13), this means dropping quantum effects associated with high energy photon emission. 
In this limit, the t-integral in (4.14) can be performed and yields a delta function, 

OO (/>2 

/*cos(t/^)=.^W.-A), (4.15) 

which says that the coherent lightfront time integrals in (4.14) become incoherent in the 
classical limit. This means that interference terms, i.e. quantum effects, drop out [13, 20]. 
(For a review of decoherence and the classical limit, see [91].) Performing one of the (/)- 
integrals leaves us with 

lim(4.10) = -r4 / drl r^ / d(f - — r^ + m^^ — ^ = ^ / « ^ , (4-16) 

which we recognise from the classical theory. (To obtain the final expression use Tvr = 
^((r — tt)^ + 771^) and change variables to u^ = (r — tt)^. The u^ integrals are elementary.) 
The extension of the above calculation to the remaining terms, (4.9) and (4.12), is direct. 
We find the following classical limits for our various terms: 

lim(free)^ = 7r^((/)) 

4> 
lim (loop + emission)^ = | Att^I^) + ^Ja'^(^7r^- ^^fc^) , (4-17) 
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lim(operator(^))^ = -An'^icf)) . 



Summing these three contributions, one finds the final result 



lim (P^) = g^((?^) as in LAD, LL, EFO . (4.18) 
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We have obtained classical RR directly from the h ^ limit of QED. The classical momen- 
tum of a radiating particle agrees with that predicted by the LAD, LL and EFO equations. 
We have made no approximation up to this point except that we work to order e^ in the 
coupling; the extension to higher orders is discussed at the end of this paper. We close this 
section with some comments on our result. 

1) The h ^ limit depends on which quantities (e.g. coupling, momenta) are chosen 
to scale with h [92, 93]. Since we wish to compare with the classical theory, where one has 
e^ but no /i, we used the 'standard' parameterisation, so for example e^ is independent of 
h, compare [92]. It is not a problem that a = e^/(47r/i) then diverges in the classical limit, 
since photon emission always brings with it a factor h. It is the cancellation of h between 
the photon momentum and the coupling which leaves a nonzero classical result. 

2) In the context of deriving classical RR from QFT, it was noted in [94] that loop 
terms could be relevant in the classical limit, see also [93]. We saw above that (loop)^ 
removes the tt^ term from p'^. Had this term not been removed, it would have caused not 
only an IR divergence, but would have been proportional to l//i, which blows up in the 
classical limit [13]. (The UV divergences we encountered were also 1/h terms). 

3) The TTjj, divergence we encountered is analogous to the classical divergence found in 
the derivation of LAD [6], which is also proportional to the momentum, or four- velocity. 
While the classical divergence is removed by mass renormalisation, it is removed by operator 
renormalisation in the quantum theory. The A:^ divergence and mass renormalisation in the 
quantum calculation are not seen classically, but are typical of lightfront quantisation. 

4) For infinite lightfront times asymptotic switching kills (operator ^■'^^) and the bound- 
ary terms, including the divergences, when integrating by parts. See [81-83] for the explicit 
calculation with exponential damping factors in nonlinear Compton scattering. (This reg- 
ularisation is gauge invariant [15].) The remaining term (4.10), from the sum of nonlinear 
Compton scattering and scattering-without-emission in Fig. 3, becomes 



{pp 



p- ■ "^^ 1 



e^ fdl_ f d^/^ to' /. , \ f . , . , '^ I ^P' ^ ^ f 7T27TI \ ,, ^, 



which agrees with the QED S^-matrix result of [13] if one neglects spin effects. (Set g = 1/2 
in [13] to go from spinor to scalar QED.) We remark that if one knew only this asymptotic 
result, and assumed it to be valid at finite time just by changing the ^-integral limits, one 
would miss (operator^-*^)) and the boundary term (4.9). These are precisely the terms which 
become Att' when h ^ and distinguish between the different classical equations^^. 

4.2 Current and position 



The above calculation can be extended to the current. The classical current for a particle 

/dr^ 
dTS\x-x,,{T))^. (4.20) 



with orbit x^^ is 



°This seems to be the reason [95] claims that S follows from QED. 
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Changing variables from r to and integrating over the spatial coordinates gives 



^ 1 f 

f = 2*;:/ **•*'• <"!' 



d(j) 

the right hand side of which can be compared to the classical limit of the current operator's 
expectation value. Here we need the full current (compare with (4.21)) 

Jm = ^\id^ - 2a^ - 2eA^)^ , (4.22) 

and A^ in this equation is as in (3.3). The calculation proceeds as above, and there is 
a ^-independent divergence which is removed by mass renormalisation. Rather than go 
through this, we present a calculation which yields the same result, but which is based on 
the position operator. As above, we calculate this operator's expectation value and take 
/i ^ 0. Localisation [96] and position [97] in relativistic QFT are difficult subjects with a 
long history, see [98], and here we are primarily interested in the classical limit. The idea 
of [55] is that, at least to lowest order, one can use the charge density Jo to define a position 
operator in instant-form quantisation by 



-/ 



d^xxJo . (4.23) 

For one-electron states with wavepackets ^2 and gi it is straightforward to show in instant- 
form quantisation that (setting x^ = for simplicity) 

f d^p / d p-^' \ 

^^^'^''^^^ ^i (2^0^^^2*^Pa-^9^ -^2^J^^^P^ ^ fe,QNw5i) , (4.24) 

in which we recognise the Newton- Wigner position operator Qnw [97]. The position op- 
erator (4.23) was applied to RR, for a different background than ours, in the papers [55, 
94, 99, 100]. Here we extend those calculations to our background field, and to lightfront 
quantisation, in which the charge density is J_ = j_. We therefore define the transverse 
and longitudinal position operators by 

Xi<P) = Jdxxj_{<P). (4.25) 

We begin with a check on the definition of the position operator by calculating the simplest 
component, obtained by replacing x in (4.25) with x+. In other words, we wish to check 
that our position operator correctly generalises to the lightfront time component. In this 
case, the contributing terms of the 'position' operator (4.25) are 



X 



Jdx j_{x+) = x+ jdp h\p)h{p) , (4.26) 



which counts the number of electrons in the system and then multiphes by x"*". Since we 
have only one electron, and net fermion number is conserved, the expectation value becomes 

{X+) =x+{ip;x+\ip;x+) =x+(in|in) = x+ , (4.27) 
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as it should be. We turn to the spatial co-ordinates, for which the expectation value is 

x~^ x~^ x^ y 

(X) = (in|X|in) + (in| f dyVi{y)X f dzVi{z)\m) - 2Re{m\X f dy f dzVi{y)Vi{z)\m) 

= : (free) + (emission) + (loop) . (4.28) 

In A+ = gauge there is no (operator^-*^^) contribution. As above, (free) is the 'free' theory 
contribution. We evaluate first the expectation value of the current in the initial state, 

(in|j_(x+)|in) =y"dpdq g*{q)g{p){p + q).ip;ix+)ip,ix+) , (4.29) 

and then perform the x-integrals in (4.25) to obtain 



(free) 



in IXI inl 



JdoWpfJ 



S-'«'^''l«'^'- 



(4.30) 



The final term gives the initial position about which the incoming wavepacket is centred at 
x+ = 0. In the limit that g also becomes sharply peaked around momentum p we obtain 



(free) = Xinit + / 



X 
kp 



A(x, = ±x(,) + 0(e^). 



(4.31) 



which are the Lorentz force results for the position and momentum of a particle in a plane 
wave, as a function of light front time (f). (As for the momentum calculation, the 'free' result 
is independent of h] the expectation value is equal to the classical result.) 

(Loop) is again a finite time one-loop self-energy contribution. Using (4.4) once more. 



(loop)-2Re^|d02|d0i|dl' 



7T2L7TI 

(kp — kl) 



Xinit + 



J kp 



4>2 
i_ r Itv 
h J kp — kl 

e ^1 



(4.32) 



The final term, (emission), comes from nonlinear Compton scattering at tree level. In the 
sequel, we present the detailed calculation of the transverse co-ordinates, setting Xinit = 
in order to to simplify the presentation. The longitudinal components may be calculated 
similarly, and we state corresponding results below. Using (4.4) again, one arrives after a 
straightforward calculation at 



(X)- 



2h^kp 



(p 2 



7T2L7T1 

kp — kl 



2 (p (p 

hj kp-kl)[j ^ J 



kl 

kp 



7T 



In this case there are no contributions from boundary terms, and no UV-divergences. (See 
[77, §11] and [101] for renormalisation of the current in QED). The classical limit is obtained 
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Figure 4. All possible cuts of these diagrams generate all possible terms to be calculated for 
(kP^) to order e^: dashed vertical lines correspond to a cut and the insertion of kP^ = kP^^\ 

as above. Taking a (/) derivative to simplify the presentation and compare directly with the 
classical calculations, the first and second lines of (4.33) contribute, respectively, 

^A7r^"(0) , ^A^^"(0) + ^|a^2(^-7r(0))" . (4.34) 

Summing these two, we find the transverse position given in (2.15) by a subset of the 
classical equations. Extended this calculation to the remaining longitudinal component, 
one finds in the end 



lim ^(X^) = ^x^((/)) as in LAD, LL, EFO, S . 
h^o d(p d(f) 



(4.35) 



Combining the two results (4.35) and (4.18), we see that, of the classical equations consid- 
ered above, only three, LAD, LL and EFO are consistent with QED to this order. Further 
(and as expected) q = mx even with recoil effects. 

4.3 Higher orders 

To distinguish between those equations which we have found to be consistent with QED, 
one needs to consider order e^ effects at finite time. (All considered equations agree on 
the post-pulse, asymptotic momentum to order e^, which can be obtained from the order 
e^ result (2.14) and Larmor's formula.) One option is to calculate k{P^) as a function of 
lightfront time and compare it to the classical result 



kp 



1 + 



A 



m^ 



/«'^^^(/"'0- 



+ C 



A2 



J2 



m^ 



(4.36) 



where Clad = 2, Cll = and Cefo = 1/2. 

We expect the next order calculation to be more difficult than that presented here, 
though. There are a few reasons for this. Consider k{P^) = k{P^^^). The two diagrams in 
Fig. 4 represents all terms in this expectation value, to order e^. A vertex is the three-point 
vertex in (3.4) (including the instantaneous photon interaction). One obtains the terms 
to be calculated, including lightfront time-ordering, from all possible cuts of the diagrams 
(red dashed lines) with the operator kP^^^ inserted at the cut. The part of the diagram to 
the left (right) of the cut then belongs to the bra (ket). This notation serves to reflect the 
unity of different contributions; in the second diagram, the loops given by the left and right 
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f^ 
^W^^ 







Figure 5. Order e^ contributions to RR come from all possible cuts of these diagrams. (The 
required counterterms are not shown.) Vacuum polarisation and pair production effects are confined 
to all possible cuts of the second and third diagrams in the third line, the latter containing the four- 
scalar interaction. 



cuts must be retained and treated together with the emission diagram, obtained from the 
central cut, in order to remove IR divergences. There are many more terms to calculate at 
order e^, obtained from all cuts of the diagrams in Fig. 5, with kP^^^ inserted at the cut. 

Most of the diagrams in Fig. 5 describe multi-photon emissions, along with loop cor- 
rections to those emissions. The loops cannot be dropped from the outset. Effects from 
pair production and vacuum birefringence are confined to the second diagram in the third 
line. Both are quantum processes, and so one expects that they should drop out as /i ^ 0, 
or that there is a parameter regime in which they are much less important than those in 
the photon emission diagrams [17, 18]. If one is only interested in the classical limit, it 
should be sufficient to drop the pair-birefringence diagram (and the final, instantaneous- 
scalar diagram in Fig. 5). We believe that 1/h terms cancel amongst themselves within 
the diagrams shown, i.e. within the groups of terms obtained from cutting the diagrams. 
At higher orders, one might also expect that UV-divergent terms and order e terms can 
combine, through 1/e * e = 1 cancellations, to give finite terms. 

One can also compare QED results with the transverse components of the classical 
orbit, dx^/dc/). We believe it would be simplest to do so using the current operator, as 
in (4.21) because in this case there are no divergences at order e^ and the calculation is 
technically simpler than when using x^ (as there are no momentum derivatives to calculate). 

In this paper we have focussed on the momentum of the electron. To see RR in the 
photon spectrum one also needs to consider at least an order e^ calculation (requiring 
emission of two photons [102], or more [17, 18]), since the photon spectrum contains no RR 
effects at order e^ [13]. This is why no RR effects were seen in previous calculations of the 
photon momentum from nonlinear Compton scattering [80-85]. 
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5 Conclusions 

We have derived dynamical quantum and classical radiation reaction from scalar QED. 
Calculating the expectation value of the momentum and position operators, we saw that 
radiation reaction appeared at order a, due to photon emission and self-energy effects. The 
h ^ limit of our expectation values gave us the orbit of a classical, radiating particle. 
Our results imply that, of the classical equations in Sect. 2, only LAD, LL and EFO can be 
consistent with QED. This is consistent with known first-order relations between these three 
equations [7, 8, 26], and with previous results on the derivation of RR from QED [51, 55]. 

Our calculation used the Hamiltonian formalism, the Furry picture, and a plane wave 
background. This choice made light front field theory the natural framework to adopt. 
Canonical quantisation of a particle in a plane wave is simple on the lightfront [60], see 
also [59], but it is only in recent years that progress has been made on instant-form quan- 
tisation of the same theory [57, 58]. (Recall that even the Lorentz orbit of a particle in a 
plane wave has no closed- form parameterisation in instant-form time x^.) 

Since expectation values are inclusive observables, all our expressions were infra-red 
finite. We found two ultra-violet divergences. The first was analogous to that found in 
the classical derivation of LAD, being proportional to the particle's momentum, and was 
removed by a multiplicative renormalisation. The second was particular to lightfront quan- 
tisation and was removed by the usual mass counterterm of lightfront QED. Here we again 
saw the suitability of lightfront quantisation for our problem; transverse dim reg [74] can 
be used without affecting the structure of the chosen background field. 

At order e^, our methods can distinguish between LAD, LL and EFO. We expect the 
calculation to be significantly harder than that presented here due to the multitude of terms 
to be calculated, and new divergent structures which can appear. 
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A Another classical equation 

The approach of [41] (S) is to begin with two coupled equations 

q = mfx — Qr ^ X — — q + Xr . (A.l) 

m 

in which q^ — m? but x^ ^ 1. So momentum, which is on mass-shell, is not proportional 
to velocity, which is off mass-shell. With the particular choice [41] 

3 47r m*^ 3 47r m^ 

one finds 

<i = h+~ (-ffQ + -^ihfq) , (A.3) 
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which looks hke LL, but without the first term, see Table 1. Thus, Sokolov's equations have 
the form (2.1) for the momentum, but not for the velocity. 
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